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We calculate the spin-Hall conductivity of a disordered two-dimensional Rashba electron gas 
within the self-consistent Born approximation and for arbitrary values of the electron density, 
parametrized by the ratio Ef/Eq, where Ef is the Fermi level and Eq is the spin-orbit energy. 
We confirm earlier results indicating that in the limit Ef/Eq S> 1 the vertex corrections suppress 
the spin-Hall conductivity. However, for sufficiently low electron density such that Ef < Eo, we 
find that the vertex corrections no longer cancel the contribution arising from the Fermi surface, 
and they cannot therefore suppress the spin current. This is instead achieved by contributions away 
from the Fermi surface, disregarded in earlier studies, which become large when Ef < Eq. 

PACS numbers: 72.25.-b, 72.10.-d, 72.20.Dp 



The spin-Hall effect, i.e., the generation of a spin po- 
larized current transverse to the direction of an applied 
external electric field, has recently raised considerable 
interest in view of its possible application in spintron- 
ics. Alongside the extrinsic spin-Hall cffectji generated 
by the spin-orbit (SO) coupling to impurities and defects, 
much theoretical effort has been devoted to the intrinsic 
spin-Hall effect arising from the one-particle band struc- 
ture of spin-orbit coupled systems m£ The initial claim 
that for a two-dimensional (2D) electron gas subject to 
the Rashba SO coupling the intrinsic spin-Hall conduc- 
tivity, <T S H, is a universal constant (ct s h = |e|/87rj^ where 
e is the electron charge) has been shown to be invalid 
even for an arbitrarily small concentration of impurities, 
which reduce <t s h to zeroA^LS At the same time, how- 
ever, for other models of SO coupling, like, for exam- 
ple, the three-dimensional (3D) Dresselhaus model, 9 the 
Luttinger model for valence band holes^S or generalized 
Rashba models taking into account non linear momen- 
tum dependences of the SO interactionii*i2ii& or a non- 
quadratic unperturbed band spectrum^4 ct sH has been 
found to be robust against non-magnetic impurity scat- 
terings, suggesting that the vanishing of <t s h is a peculiar 
feature of the linear Rashba model. This is also sup- 
ported by rather general arguments which do not rely on 
the specific scattering process ii^ 

Within the linear response theory, the vanishing of 
(J s h in the Rashba model has been ascribed to a can- 
cellation effect of the spin-dependent part of the ladder 
current vertex in the Born approximation for impurity 
scattering £&LiL This cancellation basically follows from 
the fact that, as long as the Fermi energy Ep is much 
larger than the spin-orbit energy Eq, the factor r™ 1 as- 
sociated with the current vertex (where r is the electron 
life-time due to impurities) is balanced by the factor r 
arising from the product of two Greens functions appear- 
ing in the current vertex kernel. However, similarly to 
what happens for other properties (e.g., the conductiv- 



ity in impure metals), such kinds of balance effects are 
usually peculiar to the assumption that Ep is the largest 
energy scale of the problem, and one may wonder if the 
cancellation mecha nism based on the vertex function de- 
scribed in Refs.[4ll5l l7l8| is still valid when Ep is compa- 
rable with Eq ■ The clarification of this issue is important 
not only to assess the role of vertex corrections in a gen- 
eral context, but it is also quite crucial in view of the re- 
cent progress made in fabricating systems with large SO 
couplings for which the Ep 3> Eq approximation may 
not be appropriate. Examples of such large SO systems 
are, among others, HgTe quantum wells ^ the surface 
states of metals and seminietalsj 1 ii& and even the heavy 
Fermion superconductor CePtaSiJ^ However the most 
striking example is provided by bismuth/silver(lll) sur- 
face alloys displaying quadratic unperturbed bands split 
by a Rashba energy of about Eq = 220 meVm& In this 
system the Fermi energy can be tuned by doping with 
lead atoms in such a way that Ep may be larger or lower 
than E &2L 

In this paper we investigate the spin-Hall conductivity 
in the Born approximation of impurity scattering and for 
arbitrary values of Ep/Eq. We find that, apart from the 
high density limit Ep/Eq — > oo, the spin-dependent part 
of the vertex function is generally not zero, and increases 
as Ep /Eq decreases, eventually reaching unity as Ep — » 
0. In this situation, the spin-Hall conductivity ct s h w ould 
be nonzero if calculated along the lines of Refs. |4l5B 7llaL 
in contradiction with the general arguments of Refs.0, 
n~5| . We resolve this inconsistency by showing that <t s h 
is actually canceled by the contributions away from the 
Fermi surface, which have magnitude equal to those on 
the Fermi surface, but opposite in sign. 

We consider a 2D Rashba electron gas whose Hamilto- 
nian is 

H 2 k 2 

H = — \-j(k x <7 y - k y a x ), (1) 
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where m is the electron mass, k is the electron wave- 
number, and 7 is the SO coupling. The correspond- 
ing electron dispersion consists of two branches Ef. = 
h 2 (k + sko) 2 /2m where s = ±1 is the helicity number 
and fco = rwy/h 2 . In the following, we parametrize the 
SO interaction by the Rashba energy Eq = h 2 k^/2m 
that, for a clean system, corresponds the the minimum 
inter-band excitation energy for an electron at the bot- 
tom of the lower (s = —1) band. For simplicity, we 
consider a short-ranged impurity potential of the form 
V(r) = V[ mp J2i <K r ~ where the summation is per- 
formed over random positions Ri of the impurity scatter- 
ers with density n%. The corresponding electron Green's 
function is a 2 x 2 matrix in the spin space, 

G(k,uo n ) = ~ + s(k x cr y - ky<j x )]G s (k,iuj n ), (2) 

1 s=± 

where G s {k, ico n ) = iui n — E% + \i — H(iu) n ) is the elec- 
tron propagator in the helicity basis, \i is the chemical 
potential, tu n — (2n + l)7rT is the Matsubara frequency 
where T is the temperature, and Y,{iu) n ) is the impurity 
self-energy in the self-consistent Born approximation, 



E(iw„) = 



1 



dk 



: G(k,iw n ) 
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where r _1 = 2nriiV^ a N /h is the scattering rate for a 
2D electron gas with zero SO interaction and density of 
states No — m/2irh 2 per spin direction. The spin-Hall 
conductivity is obtained from 



cr sH = - lim 



lxsxK ac (u> + iS) 



(4) 



where K sc is the spin-current-charge-current correlation 
function given by 



K sc (iu) rn ) 



r 



dk 



Tr[^(k)G(k,^„ 



Jc(k,iu n + iuj m ,iu n )G(k, iui n )} . 



(5) 



Here j%(k) = {S z ,v y (k)}/2 — h 2 k y a z /2m is the current 
operator in the y direction for spins polarized along z 
and J® is the vertex function for charge current along 
the x direction. In the Born approximation for impurity 
scattering, J* satisfies the following ladder equation: 

Jc(k,iui,iuj n ) = j*(k) 
1 f dk' 

G(k', iuj t )J*(k', iuj h iuj n )G{k',iuj n ), 

(6) 
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where Jc(k) = ev x (k) = ehk x /m + cya y /H is the 
bare charge current. Equation 10 can be rewritten as 



J* (k, iwi, iuj n ) — ehk x /m + ejT(iu!i, icu n )/h where F rep- 
resents the SO corrections to the charge current func- 
tion. By using Eq. (J2Jl and by taking advantage of the 
momentum independence of S, the correlation function 
K sc reduces to 

K sc (iuj m ) = i^—^-T^2r y (iu n +iuj m ,iuj n ) 

n 

xB^iun + iu} m ,iuj n ) 
.eh 2 j 
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T ^2 IC(iuj n + iu) m , iuj n ) (7) 



where T y is the component of T proportional to cr y , T y 
Tr(CT !/ r)/2, which, by using Eq.@), becomes 



T y (iuji,iuj n ) 



8tttNq + T-B 2 (ioJi, iuj n ) 



(8) 



8nrN - B 3 (iuji,iLu n ) 
In Eqs.JJJ and JSJ) the function B\, B2, and B3 are 

/dk 
— k 2 ^2sG- s (k,iu;i)G s (k,iuj n ), (9) 

S 

/dk , 
— A; 2 J^s G s {k, iui)G s {k, iv n ), (10) 
s 

/dk 
— k^2G s (k,iuji)G s '(k,iuj n ). (11) 

At this point, the analytical continuation to the real 
axis, iu) m — > to + id, can be performed by following the 
usual steps^ leading to 

T ^2 fcii^n + iv m , iuj n ) 



^[/(e + w) - /(e)]Im£(e + u + i5,e- i5) 

} 27T 

^ [/(e + u) + /(e)]Im/C(e + lu + id, e + id), (12) 
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where f(x) = l/[exp(x/T) + 1] is the Fermi distribution 
function. When the spin-Hall conductivity is evaluated 
via Eq.||lJ, it is clear that the resulting er sH will be given 
by the sum of two contributions, and cr^, respec- 
tively, defined as the first and second line in the right- 
hand side of Eq. , and characterized by different com- 
binations of retarded (R) and advanced (A) Green's func- 
tions (see below). The first term, cr^f, contains in the 
limit lu — > the term df(e)/de which, for T = 0, restricts 
all quasiparticle contributions to the Fermi surface. The 
second term instead has an integral containing /(e), al- 
lowing therefor e fo r processes away from the Fermi sur- 
face. In Refs.lJyjTjOl this term has been disregarded he- 
cause in the large Ep limit it scales as Eq/(tE f )^^ and 
the spin-Hall conductivity has been approximated by the 
a^f contribution alone, which at zero temperature re- 
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duces to 



RA 



eh 2 -y 



\RA 



(13) 

Here G R< " A ^ is the retarded (advanced) Green's function 
and T RA = T y (0 + iS,0 — iS) is the ladder vertex func- 
tion Jsj calculated at iu>i = iS and iui n — —iS. By as- 
suming that the SO energy E is negligible with respect 
to Ep, then the self-energy can be approximated as 
T, R (u>) = — i/2t^ and the bubble term B 2 defined in 
Ea. lfTnjl reduces to B 2 (i5,—iS) = ~8irTN k , which by 
using Eq.JSJ, leads to T RA = 0. This is the ver tex can- 
cellation mechanism pointed out in Refs.|4ll5l l7l8j ]. 

We reexamine now Ea. l|13(l by relaxing the hypothe- 
sis Ep 3> Eq. For practical purposes, we introduce an 
upper momentum cut-off k c such that all the relevant 
momentum and energy scales are much smaller than the 
corresponding cut-off quantities, namely k ,kp <C k c , 
E F ,E Q , 1/t <C E c — h 2 k 2 /2m. After the analytical con- 
tinuation, the integration over momenta in Eq.@ can 
be performed analytically and the real axis self-energy is 
evaluated numerically by iteration. The obtained £ is 
then substituted into T RA and Ea.fHfy. whose mo- 
mentum integration allows for an analytical evaluation 
due to the momentum independence of E. To explore 
the effect of varying Ep on the spin-Hall conductivity, we 
have first evaluated the Green's functions at fixed num- 
ber electron density n, where n = 2 (n = 0) means that 
all states below the cut-off energy E c are filled (empty), 
and subsequently the corresponding Ep for a given n has 
been extracted from 



1 f 00 

55 



N a (u) 
N ' 



(14) 



where N s (uj) = -(1/tt) / dk/2nkIvaGf(k,uj) is the den- 
sity of states for the interacting system and f(u>) = 
9{—uj) at zero temperature. 

In Fig. LHa) we report the SO vertex function T RA as 
a function of Eft and for several values of the SO energy 
Eq ranging from E t = 0.8 up to E t = 4 (from bottom 
to top) . The coupling to the impurity potential has been 
set equal to E c t = 80 in all cases. The corresponding 
values of the number electron density n as a function of 
Ep of the interacting system are plotted in the inset of 
Fig. ma). For Epr ~ 10, the Fermi energy Ep is suffi- 
ciently large compared to E and T RA is negligibly small, 



confirming the results reported in Refs. 4 5 7 8]. How- 
ever, as Epr is decreased, F RA increases monotonically 
up to T RA ~ 1 for Ep/Eq ~ 0. In these circumstances, 
therefore, the vertex cancellation mechanism is no longer 
active, and the corresponding spin-Hall conductivity a R ^f 
is expected to be non-zero. This is indeed shown in Fig. 
^b) where cr^j 4 , Ea. ljT3|) . is plotted in units of |e|/87r as 
a function of Epr. The nonmonotonic behavior of 
is due to the competition between the increase of T RA 
shown in Fig. Qfa) and the decrease of the integral ap- 
pearing in Ea. i|13fl as Ep —> 0. 
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FIG. 1: (a) Spin-orbit vertex function T^ A as a function of 
Eft for different values of the spin-orbit energy Eq. The scat- 
tering time t has been set equal to r = 80/ E c , where E c is the 
upper energy cut-off (see text). Inset: number electron den- 
sity n of the interacting system as a function of Eft. (b) The 
retarded-advanced part cr^j 4 of the spin-Hall conductivity in 
units of |e|/87r obtained from Eg. 1131 for the same parameter 
values of (a). 



A non-vanishing spin-Hall conductivity in an impure 
2D Rashba electrongas is at odds with the general ar- 
guments of Refs. 0113 where ct s h has been shown to be 
zero for any spin-conserving momentum scattering, in- 
dependently of the ratio Eo/Ep. However, as already 
pointed out above, the physical spin-Hall response is not 
entirely defined by Eq.lJTSJ, but should also include the 
contributions away from the Fermi surface given by the 
second term in the right-hand side of Eq. l|12|) . Hence 
c s h = o-^j 4 + a™, where 



r RR 



'sH 



eh 2 -f 



Im 



def(e)r™(e) 
dG R s (k,e 



de 



G R (k,e), (15) 



and T RR (e) = T y {e + iS,e + iS). 

Our numerical calculations of <t^j , Eci. (|15fl . are plot- 
ted in Fig. |21 (dashed lines) together with the correspond- 
ing a^f results already plotted in Fig. LT^b). For all 
Ep/Eq values, has the same magnitude of a^f but 
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FIG. 2: The different cotributions to the spin-Hall conduc- 
tivity for the same parameters of FigQ (solid lines), 
a 5P (dashed lines), and the physical spin-Hall conductivity 
<JsH = c^h 4 + clM i (gray lines). All conductivities are given in 
units of |e|/87r. 

with opposite sign, so that the resulting physical spin- 
Hall conductivity, cr S H = cr^f + (gray lines) reduces 
to zero within the accuracy of our numerical calculations. 

The results plotted in Fig[5] clearly demonstrate that, 
generally, a correct evaluation of the spin-Hall conductiv- 
ity must take into account the contributions l|15H away 
from the Fermi surface, resolving therefore the concerns 
expressed in Ref . [llj about an only-on- Fermi-surface can- 
cellation mechanism. However, on this point, a few re- 
marks should be brought to attention. First, the cancel- 
lation between a^f and a^ 1 suggests that, by suitable 
mathematical transformations, the (nominal) off-Fermi 
surface contribution (|15|l may be expressed as — cr^, re- 
sulting in a cancellation mechanism that is, after all, a 
Fermi surface property. However, we have been unable 
to find such a transformation. A second possibility is 
that <7^p is, generally, a genuine off-Fermi surface quan- 
tity, but that, accidentally, for the model hamiltonian of 
Eq. JJJ, such a term is quantitatively equal to — cr^. In 
this case, any variation from the linear Rashba model of 
QJ would result in and a^p terms which do not 



mutually cancel, leading to a nonzero spin-Hall conduc- 
tivity. In this respect, one should note that, in fact, the 
general arguments of Refs.00] about the vanishing of 
ct s h apply only for model hamiltonians of the type l|T[l. 

Before concluding, it is worth stressing that the results 
presented in this work could be relevant also for systems 
described by non-linear Rashba or 3D Dresselhaus SO 
couplings, or by non-quadratic unperturbed electronic 
band structures. It is known that for such systems, the 
spin-Hall conductivity in the presence of momentum scat- 
tering is non-zero also for Ef/Eq — > oo because the SO 
vertex does not vanish AiSiiiiiSiiii^ Our results suggest, 
however, that even in this case, for finite Ep /Eq, a quan- 
titatively reliable calculation of ct s h should take into ac- 
count also the off- Fermi surface contributions. This could 
be for example the case of the system studied in Ref. |isj 
where Ep is of the same order as Eq and the unperturbed 
band spectrum is clearly non-quadratic. 

In conclusion, we have calculated the spin-Hall conduc- 
tivity <7 s h for a 2D electron gas subjected to the linear 
Rashba SO coupling in the Born approximation for im- 
purity scattering. We have shown that, apart from the 
Ep — ► oo limit, the spin-dependent part of the vertex 
function is nonzero and increases as Ep — > 0, leading to 
nonzero Fermi surface contribution to the spin-Hall 
conductivity. We have demonstrated that the physical 
spin-Hall conductivity ct s h actually includes also contri- 
butions away from the Fermi surface, erjfp, which are as 
large as those on the Fermi surface, but of opposite sign, 
leading to a vanishing ct s h for arbitrary values of Ep / Eq . 
We expect that, given the arguments of Refs. |7lllq . the 
mutual cancellation of er^ 4 and for Ep < oo holds 
true also beyond the self-consistent Born approximation 
employed here. 
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